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Abstract. We report on observations of correlated behavior between the prompt y-ray and optical emission from GRB 
0803 19B, which confirm that (i) they occurred within the same astrophysical source region and (ii) their respective radiation 
mechanisms were dynamically coupled. Our results, based upon a new cross-correlation function (CCF) methodology for 
determining the time-resolved spectral lag, are summarized as follows. First, the evolution in the arrival offset of prompt 
y-ray photon counts between Swift-BAT 15-25 keV and 50-100 keV energy bands (intrinsic y-ray spectral lag) appears to be 
anti-correlated with the arrival offset between prompt 15-350 keV y-rays and the optical emission observed by TORTORA 
(extrinsic optical/y-ray spectral lag), thus effectively partitioning the burst into two main episodes at ~ T + 28±2 sec. 
Second, the rise and decline of prompt optical emission at ~ T + 10 ± 1 sec and ~ T + 50 ± 1 sec, respectively, both coincide 
with discontinuities in the hard to soft evolution of the photon index for a power law fit to 15-150 keV Swift-BAT data 
at ~ T + 8 ± 2 sec and ~ T + 48 ± 1 sec. These spectral energy changes also coincide with intervals whose time-resolved 
spectral lag values are consistent with zero, at ~ T + 12 ± 2 sec and ~ T + 50 ± 2 sec. These results, which are robust across 
heuristic permutations of Swift-BAT energy channels and varying temporal bin resolution, have also been corroborated via 
independent analysis of Konus-Wind data. This potential discovery may provide the first observational evidence for an implicit 
connection between spectral lags and GRB emission mechanisms in the context of canonical fireball phenomenology. Future 
work includes exploring a subset of bursts with prompt optical emission to probe the unique or ubiquitous nature of this result. 
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1. INTRODUCTION 

Swift's unique dynamic response and spatial localization precision, in conjunction with correlative ground-based 
follow-up efforts, has resulted in the collaborative broad-band observations of GRB 0803 19B [1]. In the context 
of an analysis focused on confronting the lag-luminosity relation yj] in the Swift era via the Burst Alert Telescope 
(BAT) i], a correlation was observed between the evolution of time-resolved spectral lag and the behavior of the 
extraordinarily well-sampled prompt optical emission light curve associated with GRB 0803 19B [4]. In general, the 
spectral lag is determined via either a peak pulse fit [5] or cross-correlation function (CCF) analysis [$]. Previous 
studies have reported on the variability of spectral lag throughout burst emission J7D, as well as its correlation to pulse 
evolution [8]. In this work, we develop a new method to calculate the time-resolved spectral lag of GRB 080319B, 
within 10 emission periods, as illustrated in Figure [T] via a modification to the traditional CCF approach. In this 
manner, we are able to explore for the first time the evolution and correlation of the intrinsic y-ray spectral lag with 
prompt optical emission. We interpret these correlated behaviors as strong observational evidence that the prompt 
optical and y-ray emission of GRB 0803 19B took place within the same astrophysical source region, with indications 
that their respective radiation mechanisms were dynamically coupled throughout the prompt phase of the burst. 
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FIGURE 1. 64-ms binned light curves for Swift-BAT 15-25 keV (Panel A), 25-50 keV (Panel B), 50-100 keV (Panel C), 100-300 
keV (Panel D) and Konus-Wind 300-1160 keV (Panel E) energy band passes. Gray regions indicate the sample periods for time 
resolved analysis. BAT error bars (Panels A-D) indicate lo" uncertainty. Konus-Wind data was taken from a supplement online Q)]. 



2. METHODOLOGY 

The Burst Advocate (BA) script was used to generate BAT event data and quality maps, which were used to construct 
32 ms binnec0 (mask-weightecQ, background-subtracted) light curves (temporal photon spectra) for = 15 — 25 keV 
and Eb = 50— 100 keV energy band^J via the BATBINEVT analysis task. The resultant F I T s files were analyzed via 
a modified CCF method [6] in order to quantify the temporal correlation between the two series of GRB light curves 
in differing energy channels, i.e. Ea and Eg, via the following: 

CCF( to , VEA , VEa ) = ^ V ^ +t »K cCF AB , (!) 

<*VE A a VE B 



3 A binning of 32 ms was used based upon a heuristic exploration of various bin resolutions (e.g. 4 ms, 8 ms, 16 ms, 32 ms, 64 ms, etc.). Once a 
spectral lag measurement was made, it was reproduced using light curves of several adjacent bin resolutions. In principle, the binning of the light 
curve should be smaller than the spectral lag in order to resolve the CCF peak. 

4 Although one increases the signal-to-noise ratio (SNR) using non-mask-tagged (raw, non-background subtracted) light curves, artifacts either 
intrinsic to the detector or associated with the background may mimic signal pulses and thus be treated as such in a CCF analysis. Also, the 
background may not be trivial to model since it may be highly variable due to slewing. Hence, in order to be conservative and consistent with future 
studies, background subtracted, i.e. mask-tagged, light curves were used throughout the analysis. 

5 Permutations of canonical BAT channel (1-4) light curve pairings were investigated. Ultimately BAT channels 1 (15-25 keV) and 3 (50-100 keV) 
were used in the analysis since they represented the largest energy differential with the highest count rate significance. 
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FIGURE 2. Swift-BAT prompt y-ray (32 ms bin) light curves for canonical energy channels 1 (15-25 keV) and 3 (50-100 keV), 
with corresponding CCF as a function of spectral lag time. The peak of a Gaussian fit (solid line) determined the time-resolved 
spectral lag. The results for sample time intervals 4, 6 and 7 (defined in FigureQJ, include 6±4 ms (Panel A), 65 ± 15 ms (Panel B) 
and 149 ± 57 ms (Panel D), respectively. Note that the time-resolved spectral lag of a given interval (Panels A, B & D) can differ 
greatly from the time-averaged spectral lag of 29 ± 4 ms, as illustrated in Panel C. All error bars denote 1 a uncertainty. 



where ov £ „ = ( v e„) 1 ^ 2 - The CCF between and Eg peaks at a given temporal offset (f„), known as the spectral lag 
(tab), which is defined as positive if the systematic shift in the arrival times of photon counts between pairs of light 
curves results in higher energy (Eb) photons arriving before those of lower energy (Ea). The generic CCF is based 
upon the Pearson Correlation (Equation [TJ. If the mean is subtractecfl then v(t) = d(t) — (d(t)). If the mean is not 
subtracted, then v(f ) = d(t) — b(t), where an attempt is made to remove background counts b(t). Our tests indicated 
that both methods agreed when one calculates the time-averaged spectral lag, i.e. the spectral lag determined over the 
entire duration of the GRB (Figure|2] Panel C). This included a test of subtracting the mean of a light curve that was 
nested within two large background intervals, which had the effect of reducing the mean and increasing the correlation 
amplitude. Note that if one takes an infinite interval, the mean would go to background, which fluctuates about zero, 
resulting in effectively not subtracting the mean in the first place, as prescribed for transient sources such as GRBs @]. 

However, when one tries to extract a time-resolved spectral lag, i.e. the spectral lag over a segment of the light curve, 
then only the mean subtraction method consistently works. This is a consequence of edge effects introduced by the 
assumptions of a stationary versus a transient temporal signal characterization, which destroy/diminish the intrinsic 
lag. If one subtracts the mean, there is an assumption that the signal exists as a stationary wave throughout time, i.e. 
beyond the sampled time series. If one does not subtract the mean, then the assumption is that the signal is entirely 
contained within the time series. Both methods disregard data beyond the overlapped interval displaced via bin shifts. 

Hence, for a bright GRB like 0803 19B, a non-mean subtracted time-resolved spectral lag analysis is dominated by 
the segment of the light curve described by a square wave, and the real light curve signal is treated as fluctuations. A 
mean-subtracted time-resolved spectral lag analysis will extract the spectral lag since it smoothes out the amplitude 



6 Mean subtraction is the formalism identical to that used by the C_CORRELATE function of IDL 6.2. The formal definition may be found online 

(http : //cl . dmf . arm. gov/base/idl_6 . 2 /C_CORRELATE . html). 



and focuses on the signal variance. Hence, we utilized the mean-subtracted definition of the CCF, since it is both 
more conservative and robust. Although it has been demonstrated that spectral lag variability is ubiquitous in GRBs 
10, [U, previous methods required relatively well-behaved sub pulses, such as those canonically described as FREDs 
(fast-rise, exponential decay). Our new methodology treats structure without regard to its functional form and affords 
greater flexibility for time-resolved spectral lag evolution studies, as illustrated in Figure [2] which contrasts sample 
periods 4 (Panel A), 6 (Panel B) and 7 (Panel D) with the time-averaged interval (Panel C) for GRB 0803 19B. 

In order to test the viability of an interval, a large known spectral lag was artificially introduced, which shifted 
one of the test light curve pairs. Intervals where the analysis was unable to recover the artificial spectral lag due to 
either low SNR or lack of structure were not used. The CCF method, regardless of definition, is sensitive to the SNR 
and, more importantly, to the shape of the region. Morphological tests using simple square, triangular and Gaussian 
shaped pulses have illustrated that both mean and non-mean subtracted CCF definitions failed to recover the artificial 
spectral lag that was obscured in regions where (i) only slopes were compared, (ii) SNR was low and/or (iii) error bars 
were large. This puts a fundamental limit on determining the time-resolved spectral lag for discrete intervals, which 
precluded an arbitrary systematic sampling of GRB 0803 19B, as illustrated in Figures[TJand [3](Panels C & D). Hence, 
although our methodology is robust, it does require a selection on intervals with structured variability. This is not an 
issue for time-averaged results, since the burst is nested between periods of quiescence, as given in Figure [2] Panel C. 

Considerable effort was given to ascribing a statistical significance to the spectral lag evolution via the determination 
of errors. Several methods were explored which included an approximation of confidence intervals, discrete CCF 
formulations and parameterizations accounting for the total height (h), half-width at half-maximum (W e ) and number 
of light curve bins (n) for a given CCF peak H], such that: ff T « i+l'fa^i) 1 / 2 ' 

Ultimately, we implemented a more intuitive approach, which featured a double Monte Carlo (MC) schema via 
canonical bootstrapping techniques fioll . In this manner, data-based simulations were used to estimate the standard 
error via the standard deviation of the bootstrapped light curve replications. For a given light curve bin {LCtin ± &LC bin ) > 
a simulated realization [LC' bin ± <7tc fc;n ) of the count rate was constructed as follows: LC' hin — LCtin + <= c7 lc w „- The 
random number seed (<!;) was determined via the IDL RANDOMInQ function, which is based upon the Box-Muller 
method and returns a normally-distributed, floating-point, pseudo-random number with a mean of zero and a standard 
deviation of one. 

The error bars on the CCF function were determined via a primary MC simulation, which were constructed from 
1000 realizations of light curve pairs, within the la error bars of the data. The distribution of CCF values within a 
given bin defined its error bar, as illustrated in Figure [2] Once the error bars were fixed in this manner, a secondary 
MC simulation constructed a subsequent set of 1000 realizations of light curve pairs, within the Iff error bars of the 
data. Each resulting CCF (with errors fixed from the primary MC simulation) was subjected to a Gaussian fi{^| 

y=Ae- {x - x < )2/ ( 2w2 ), (2) 

whose peak was extracted via the IDL function MPFITPEAfQ. The peak fit to the data was used to determine the 
spectral lag value, while the numerical standard deviation of the distribution of simulated peaks was used to determine 
the error, which was in agreement with comparisons made with the cumulative fraction at 68% and CCF shape (ff T ), 
to within a factor of ~ 0.80 ± 0.26. We interpret this as the 1 ff spectral lag erroiFl as illustrated in Figure[3] Panel D. 

The primary MC simulation optimized the CCF value per bin and provided a straight forward generation of CCF 
error bars, which are typically omitted, without appealing to cryptic statistical arguments or the assumptions of 
otherwise un-weighted fits. The secondary MC simulation optimized the CCF peak and simulated our data analysis. 
Since the CCF is normalized to unity, the correlation coefficient ascribes the confidence level of the correlation, which 
ranged from ^40%-90%, with typical values of ^60%. Segments with correlation values <~ 30% were omitted. 

In general, spectral lag calculations are sensitive to a series of selection parameters that include the energy band pass 
of each comparative light curve, temporal bin resolution and emission interval. The CCF and pulse fit methods are 
both based upon a functional expression for the peak fit (Gaussian, cubic, etc.), while the former requires a bin shift 



7 RANDOMN is described online at http://idl astro. gsfc.nasa. gov/idl_html_help/RANDOMN . html. 

8 Typically the CCF is asymmetric, with a distribution typically skewed to the right. This is due in part to the intrinsic hard to soft pulse width 
evolution. The CCF width is also affected by the smearing of multi-pulses within the segment. The skewness was not commonly minimized via the 
introduction of error bars from the primary MC simulation. In such cases, a region of local symmetry about the peak was selected using asymmetric 
endpoints. For extreme cases, a cubic function may be used to guide the functional fit in order to accommodate the skewness. 

9 See http : / /purl . com/ net /mpf it, for more information on MPFIT. 

10 We note that taking the peak fit error to the data resulted in values that were smaller by a factor of ~ 2.4 ± 0.80. 
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FIGURE 3. Panel A - TORTORA prompt optical flux density for ~1.3 second exposure intervals (solid black circles) Qj], A 
cumulative fit (solid line) is based upon three Gaussian curves with peaks at ~T+19 sec (dashed line), ~T+34 sec (dotted line) and 
~T+45 sec (dash-dotted line). Panel B - Swift-BAT prompt y-ray light curves for (32 ms bin) canonical energy channels 1 (15-25 
keV, light gray line), 3 (50-100 keV, dark gray line) and (1.3 second bin) 15-350 keV (black line). A cumulative fit (solid line) is 
based upon three Gaussian curves with peaks at ~T+17 sec (dashed line), ~T+40 sec (dotted line) and ~T+49 sec (dash-dotted 
line). Panel C - Photon index (E~ a ) for power law fit to 15-150 keV Swift-BAT data, with lcr error bars. Panel D - Time-resolved, 
intrinsic spectral lag between Swift-BAT 32 ms bin light curves for 15-25 keV and 50-100 keV energy bands (solid black circles) 
and Konus-Wind 64 ms bin light curves for 18-70 keV and 70-300 keV energy bands (open gray circles). 



range. A self-consistent comparison, which includes uncertainties, must account for these selection effects. Efforts 
to standardize the spectral lag analysis are ongoing and would facilitate direct comparative studies across detector 
sensitivities and energy band passes |3fl. 



3. RESULTS & DISCUSSION 

Our preliminary results are illustrated in Figure[3] Panels A-D. A major result from our analysis is that the first episode 

in spectral softening of the photon index (Jjj^ £ y a j at ~ T + 8 ± 2, based upon time-resolved power law fits to 

15-150 keV Swift-BAT data intervals, coincides with both the initial rise of the optical flux and first apparent peak 
in spectral lag at ~ T + 10 ± 1 sec, as illustrated in Figure[3] Panels A, C and D. This is consistent with independent 
Konus-Wind time-resolved spectral analyses of hardness ratios [1]. Although the onset of the afterglow may obscure 
the analysis of prompt emission for t >^ T + 50 sec, the second episode of spectral softening at ~ T + 48 ± 1 
sec coincides with the apparent decline of the optical flux and zero spectral lag at ~ T + 50 ± 2 sec, which is also 
independently corroborated via Konus-Wind time-resolved analyses [1], as illustrated in Figure[3] Panels A, C and D. 



The generic agreement of the overall temporal coincidence and morphology between the prompt y-ray and optical 
light curves (Figure [31 Panels A & B) in conjunction with the observed spectral and temporal evolution (Figure [5] 
Panels C & D), suggests that they arose from a common source region (cf. UHo . in agreement with the double-jet 
hypothesis where both emissions occurred within the narrow jet (JJJ _ However, separate radiation mechanisms were 
most likely responsible since the extrapolated y-ray flux density to the optical band was deficient by ^4 orders of 
magnitude when compared to observation fH [l2l [l3ll . The steep rise/decline, short duration and lack of increasing 
pulse width of the prompt optical emission disfavor external forward/reverse shocks. Hence, internal shocks have 
been suggested as the source region with synchrotron emission responsible for the optical and inverse Compton 
scattering/synchrotron self Compton for the y-ray s, with associated GeV photon emission |Ql[l2l- Alternatively, it has 
been suggested that non-relativistic forward internal shocks generated the prompt optical emission, while relativistic 
reverse internal shocks generated the prompt y-rays, with sub-GeV/MeV photon emission fl3ll . Such high energy 
emission may be tested on future GRBs using Fermi (formally known as GLAST) via joint Swift-BAT analyses 11411 . 

Furthermore, there are several observations that lead to effectively separating the burst's duration into two main 
episodes partitioned roughly at the midpoint of ~ T + 28 ± 2 sec. The first is that the (bimodal) evolution of the 
intrinsic y-ray spectral lag increases att >^ T + 28 ±2 sec, which appears to be anti-correlated with the extrinsic 
(optical/y-ray) spectral lag, which is of the order of a few seconds when t<^T + 28±2 sec (Panels A & B), as 
observed via the smoothed Gaussian fits, as illustrated in Figure[3] Panels A & B. Beyond this common midpoint, the 
optical and y-rays do not correlate as well or at least are ambiguously correlated, i.e. extrinsic (optical/y-ray) spectral 
lag is either zero or negative at later times. The ambiguity is due to apparent peak misalignment. Hence, the intrinsic 
time-resolved y-ray spectral lag is maximum att>^T + 28±2 sec, while the extrinsic time-resolved (optical/y-ray) 
spectral lag is maximum att <^ T + 28 ±2 sec. In addition, an independent analysis of BAT 15-150 keV light curves 
has revealed that the characteristic variability timescale of GRB 0803 19B was ^100 ms for t < ~ 7+28 sec and ~1 
sec for t > ~ 7+28 sec 11511 . Furthermore, independent time-resolved spectral analysis of Konus-Wind data illustrated 
that peak energy decreased from 75 1 ± 26 keV to 537 ± 28 ke V at ~ T + 24 ± 2 sec QJJ] . 

In conclusion, we interpret these correlated behaviors as strong observational evidence that the prompt optical and y- 
ray emission took place within the same astrophysical source region with dynamically coupled radiation mechanisms, 
which until now has only been conjecture JH [TH EH . This potential discovery may provide the first observational 
evidence for an implicit connection between spectral lags and GRB emission mechanisms in the context of canonical 
fireball phenomenology. A full theoretical analysis of this result is currently in preparation. Future work includes an 
application of our methodology to observations of a subset of bursts with prompt optical emission to probe either the 
unique or ubiquitous nature of this result. Ultimately, understanding the mechanism(s) responsible for spectral lag may 
reveal a fundamental and unprecedented view from within the GRB fireball and its progenitor(s). 
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